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ABSTRACT 



Aims. Planet formation has been known for many years to be tied to the spatial distribution of gas and dust in disks 
around young stars. To constrain planet formation models, imaging observations of protoplanetary disks are required. 
Methods. In this framework, we have undertaken a mid-infrared imaging survey of Herbig Ae stars, which are pre-main 
sequence stars of intermediate mass still surrounded by a large amount of circumstellar material. The observations were 
made at a wavelength of 20.5 /im with the CAMIRAS camera mounted at the Cassegrain focus of the Canada France 
Hawaii Telescope. 

Results. We report the observations of three stars, HD135344, CQTau and HD163296. The circumstellar material 
around the three objects is spatially resolved. The extensions feature a disk like shape. The images provide direct 
information on two key parameters of the disk : its inclination and its outer radius. The outer radius is found to be 
quite different from the one deduced from disk models only constrained by fitting the Spectral Energy Distribution of 
the object. Other parameters of the disk, such as flaring, dust mass have been deduced from fitting both the observed 
extension and the spectral energy distribution with sophisticated disk models. 

Conclusions. Our results show how important imaging data are to tighten constraints on the disk model parameters. 
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1. Introduction 

The formation of circumstellar disks is a natural outcome 
of the star formation process by which a molecular core col- 
lapses to form a star ijShu et alJl!987j) . Circumstellar disks 
k> , can outlive the period during which stars form and still be 
j_J ■ present when the star is in its Pre-Main-Sequence (PMS) 
C\3 phase. In these disks, composed of gas and dust, various 
physical processes can lead to the growth of dust grains 
and eventually to the formation of planets. Understanding 
the physical conditions that prevail in these objects is of 
crucial importance when studying planet formation. 
The study of circumstellar disks is a field in fast devel- 
opment both fro m the observa tional and the modeling 
point of view (e.g. iNattal l)2004|) and references there-in). 
The Infrared Space Observatory (ISO) has given clues 
on the d ust composition of a sample of isolated HAeBe 
systems l|Bouwman et al.ll2001l fMeeus et al.ll200lD . While 
these spectra reveal the composition of the dust, no direct 
information concerning the spatial distribution of the dif- 
ferent dust species can be inferred from the ISO data. Most 
studies so far have used the Spectral Energy Distributions 
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(SEDs) to put constraints on the spatial distribution of the 
circumstellar material. 

Models of protoplanetary disks are increasingly success- 
ful at accounting for much of the observed properties. For 
instance, they can justify that disks SEDs are generally 
rath er fiat in vF v , where v is the frequency and F v the 
flux ijKenvon &: Hartmannl ll987b Furthermore, models can 
explain that dust features are almost all seen i n emission 
i|Catvet et all l|1991|blChiang fe Goldrricbl (|1997|) . hereafter 
CG97), the presence of a near-in f rared e xcess in the SEDs 
of He rbig Ae stars ijNatta et al.l l|2001|) . iDullemond et all 
(2001), hereafter DDN01) and i nterpret the differences 
observed in the far-IR exce sses l|Dullemond et all (2002; 
IDullemond fc Dominikl 12004(1 . Fitting the SED only allow 
to make conjectures on how the disks' geometry looks like; 
spatially resolved imaging data of those disks are absolutely 
necessary to verify theories and models' assumptions. For 
instance, key parameters, such as the disk surface density 
profile with radius, are still very poor constrained when fit- 
ting SEDs. 

Mid-infrared imaging observations from a large ground- 
based telescope are potentially well suited to bring spa- 
tial information on disk around Herbig Ae (HAe) stars. 
HAe stars represent the middle stage of PMS evolution 
of intermediate-mass stars (~ 2-3 Mq); they are bright 
enough to heat sub-micron dust grains at 100 AU to a tem- 
perature of about 150 K. Grains at such a temperature have 
their peak of thermal emission in the mid-InfraRed (mid- 
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IR). The diffraction limited angular resolution achievable 
with a 3.6 meters class telescope in the mid-IR, 0.6/f.2 
arcsec at 10/20 /*m, corresponds to a distance of 60/120 
AU for a star located at a typical distance of 100 pc and 
thus allows a relatively good sampling of disks whose sizes 
range in the several hundreds of AU. The first attempts 
to resolve the spatial structure of the circumstellar mate- 
rial around HAe stars were performed with multi-aperture 
observations with a single bolometer; emission extending 
up to large distances from the star (more than 5 arc- 
sec) were fo und around 3 obje cts HD97048, HD97300 and 
HD176386 ijPrusti et alJll994|) . Such extended emissions 
were attributed to the emission from large molecules and 
small grains transiently heated by star-light and distributed 
in a dust shell with a large inner radius, probably a rem- 
nant from the cloud from which the star was born. It is 
only with the advent of mid-IR cameras that the full po- 
tential of mid-IR observations to study disk structures was 
achieved. Single dish observations have revealed the disk 
structure in the 500 AU ( Astronomical Units) range around 
two H Ae stars: AB Aur (Ma rsh e t al.lll995i iPant in et alJ 
120051) and HD100546 llOradv et al.M20011 ILiu et alJl2003D . 
Interferometric observations have allowed to probe the in- 
nermost regions of t he disk (1-10 AU) and numerous objects 
have been observed iMill an-Gabet et al.l2001j:lTuthill et al 
| 2002tlWilkin fc AkesorJl2003HLeinert et alJl2004ULiu et al. 
2005J. When mid-IR single dish observations have now 
clearly demonstrated on a few examples that they can pro- 
vide unique information on the disk structure, the con- 
strains brought on the disk modeling by such spatial in- 
formation has not yet been fully exploited. 

In this paper, we were interested in both increasing the 
number of HAe for which information on the mid-IR spa- 
tial extension is available, and in using this information, in 
combination with the already existing SED measurements 
between 1 and 100 /mi, to constrain the parameters of the 
disk models developed these last years. The paper is orga- 
nized as follows : the observations and data reduction are 
described in Sect. 2. Sect. 3 deals with the results, in terms 
of spatial extensions. In Sect. 4, a first, simple approach is 
used to establish the disk inclination, and in Sect. 5, we 
describe the 2-D radiative transfer code used to reproduce 
the observations. A discussion of the results follows in Sect. 
6. Conclusions and perspective are drawn in Sect. 7. 



2. Observations and data reduction 

We have observed a sample of three Herbig Ae stars: 
HD 135344, CQ Tau and HD 163296. Table □ presents 
the main stellar parameters of the sample. The ob jects 
were selected from the catalogue of iThe et all (|1994f) and 
Mal fait et alJ l)1998f) according to the following criteria: the 
objects are bright in the mid-IR, relatively close and iso- 
lated, i.e not associated to extended diffuse emission due to 
the parental cloud. 

CQ Tau is loc ated at a distance of lOOli? pc and has 
an age of 10 Myr JNatta et al.l|200ll | van den Ancker et alJ 
Il998|) . HD 163296 is at a distance of 122+ }g pc and has 
an age of 7 ± 5 Myr ijvan Boekel et afll2005|) . The distance 
and age of HD 1353 44 are more contro versial. Until 2001, 
a distance of 84 pc jMeeus et al.ll2001[) with an age of 17 
± 3 Myr llThi et aLlboOlD were use d for this object. But, 
in a recent paper jvan Boekel et al.l2005|) . the distance was 



re-evaluated to 140±42 pc and the age to 8 ± 4 Myr; these 
latter values will be used in the following. 

The observations were performed w ith the mid-IR 
camera CAMIRAS lILagage et all Il992h installed as a 
visiting instrument at the Cassegrain focus of the Canada 
France Hawaii Telescope (CFHT). The camera is equipped 
with a Boeing 128x128 pixels Blocked Impurity Band 
(BIB) detector sensitive up to a wavelength, A, of ~ 27 
/im. A filter centered at 20.5 /im with a Full Width Half 
Maximum (FWHM) bandpass, AA, of 1.11 /im was used. 
The Pixel Field of View (PFoV) on the sky was 0.29 arcsec 
; such a PFoV provides a good sampling of the diffraction 
pattern which is of 1.5 arcsec FWHM. 
The objects were observed between 2000 March, 18 th and 
2000 March, 24 th . During the run, seeing and weather 
conditions were extremely favorable and stable in time. 
HD 163296 and CQ Tau were observed at a median airmass 
of 1.4, and HD 135344 at an airmass around 1.9, which is 
the best achievable when observing from CFHT. Standard 
chopping and nodding techniques were applied to suppress 
atmosphere and telescope background emissions; the chop- 
ping throw was 16 arcsec to the North and the frequency 
used was 3.33 Hz; the nodding amplitude was 20 arcsec 
to the West. The nodding direction was perpendicular to 
the chopping direction, in order to get the best spatial 
resolution; given the low chopping and nodding throw and 
the field of view of the camera, the source always remained 
within the detector field of view, the obtained images 
contain thus 4 beams (2 positive, 2 negative). Given the 
huge photon background in the mid-IR, the elementary 
integration time was set to 15 ms, and the images were 
co-added in real time in order to store only two co-added 
images (one for each chopping position) every second. 

The basic data reduction is standard. The data cubes 
of one observation are carefully stacked with rejection of 
corrupted planes. A shift-and-add procedure is applied to 
each cube of images using a correlation based method with 
a re-sampling factor of 8:1. The four beams are then com- 
bined in one image by a source extraction algorithm fol- 
lowed by a shift-and-add procedure. Finally, flux calibra- 
tion is achieved via aperture photometry of a set of photo- 
metric standard sta rs such as aTau, aBoo, /3Gem or 7Dra 
l|Cohen et al.lll999(l . The photometry gives a total flux of 
6.3±0.6 Jy for HD 135344, 23±3 Jy for CQ Tau and 18±4 
Jy for HD 163296, in good agreement with the IRAS values 
(Tab.HJ. 

3. Extended emission 

The three objects of the sample are spatially resolved. This 
can be seen on Fig. ^ where we have compared the average 
annular profile of the object with those of the Point Spread 
Function (PSF), obtained from the observation of a refer- 
ence point-like star. Extended emission is detected up to 
100-300 AU. 

We carefully checked that the observed extensions are 
not artifacts, but the result of true extended emission from 
the objects. Several arguments lead us to reject explana- 
tions of the extensions in terms of temporal variations of 
the PSF between the observation of the object and the ob- 
servation of the reference star. One possible cause of such 
temporal PSF variations could be variations of the seeing. 
This hypothesis is however rejected for the following rea- 



C. Doucet et al.: Mid-InfraRed imaging of Herbig Ae stars 



3 



Table 1. Properties of the three Herbig Ae stars observed with the CAMIRAS mid-InfraRed camera; first two columns: 
coordinates of the objects (right ascension and declination); third column: their spectral type; fourth column : their 
distance, as deduced from Hipparcos data; fifth column: flux in the IRAS 25 /im band; sixth column : flux at 20 
(im obtained i n this paper ; seventh column : age of the o bjects in My r s. Ref erenc es: ( 1 ) [CcmMQlL^Lal 
Dominik et all \200& ) . (3)\Javawardhana et all \200A) . (A)\Meeus et all 11200 A) . f 5)\Manninas & Saraen 



1998 



199 



(2) 
M 



Mannings & Saraeni "ttoo'd). f 7)\Natta et all \20ol). 'f8)^Svtvester et all ' \199<( ).' ' (9) tTesti et all \200\ ). (lOi TThi et all 
"200 A) . (ll) \van Boekel et all \200$) . fl2) \van den Ancker et all \1998\) . 



Object 


RA (2000) 


Dec (2000) 


Spectral type 


Distance (pc) 


F 25 (Jy) 


Fso (Jy) 


age (Myrs) 


References 


HD 135344 


15 15 48.4 


- 37 09 16 


F4V 


140+- 


6.7 


5 ± 1 


8±4 


1,2,4,8,10,11 


HD 163296 


17 56 21.4 


- 21 57 20 


A3Ve 




21 


18 ± 4 


7±5 


3,4,5,11,12 


CQ Tau 


05 35 58.4 


+24 44 54 


Al-F5IVe 


loot?? 


20.6 


23 ± 3 


10 


6,7,9,12 
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Fig. 1. Annular averages intensity aperture for the source (solid line) and PSF (dashed line) normalized to the peak value, 
as function of radius. The error drawn are 1 cr-errors RMS. The reference star for HD 135344 is a Boo, for HD 163296 
7 Dra and for CQ Tau, f3 Gem. The plot labeled '2 PSF' shows two reference stars a Boo and 7 Dra as observed on the 
24 t h of March 2000. Although the two reference stars have very different fluxes (7Dra with 43 Jy and a Boo with 170 
Jv. lVan Malderen et alJ l)2004|) '). the two profiles are similar. 



sons. First, concerning the limitations of the spatial resolu- 
tion, the seeing contribution at 20.5 /jm is negligible with 
respect to the diffraction : for a typical seeing value of 0.8" 
FWHM in the visible range, one can estimate a seeing con- 
tribution at 20.5 /im around 0.4", when using the A^ -1 / 5 -* 
scaling law; thus, seeing induces PSF changes of the order 
of 1 pixel FWHM, which is much smaller than the widths of 
observed extensions. Secondly, the seeing was quite stable 
during the observations; thus we estimated seeing variation 



effects to be much lower than one pixel. Note also that 
two of the objects were observed during different nights. 
HD 135344 was observed for four different nights: on the 
March, 18 th , 2000 (exposure time 6 mn), on the March, 
19 th (exposure time: 3 mn), on the March, 21 s * (exposure 
time: 6 mn) and on the March, 24 th (exposure time: 2 mn 
30). Its extension is confirmed over the 4 nights. CQ Tau 
was observed during 2 different nights on the 20 th of March 
and on the 21 s * of March (respectively with 3 and 6 mn of 
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integration time) and is spatially extended in both datasets. 
HD 163296 was observed only once (exposure time: 3 mn) 
on the 24 th of March. 

Another possible source of fake extensions could be 
chromatic effects. Indeed, HAe stars have large infrared 
excesses, and thus have SEDs quite different from that of 
the PSF reference stars. Any filter leak, either on the blue 
or on the red side of the nominal filter bandpass, would 
then potentially lead to PSF variations between point-like 
HAe objects and PSF reference stars. We double checked 
the filter transmission at the operating temperature of 
the filter (10K) using a Fourier Transform Spectrometer 
with a spectral resolution of 4 cm -1 . The rejection rate 
outside the filter bandpass is typically better than 10 -3 . 
We simulated the PSF variations due to such a filter 
using the mid-IR sp ectra of the objects obtained with ISO 
ijMeeus et al.ll200l|) for two of the three extended objects; 
(CQ Tau was not observed by ISO SWS). The simulated 
PSF variations lead to some extended emission, but much 
fainter than observed, both in intensity and in spatial 
extension. Thus, explanation of the observed extensions in 
terms of chromatic variations of the PSF can be discarded. 

Note also that not all the observed HAe stars observed 
are extended. In order to make comparison, a far away HAe 
star, with a simi lar (even more) I R excess, HD179218 lo- 
cated at 240 pc ijMeeus et all200i|) , was observed. This ob- 
ject does not show extended emission at 20.5 /im (Fig. 
This is an additional argument to conclude that the exten- 
sion observed in HD 135344, CQ Tau, HD 163296 is really 
due to an extended emission from the objects. 

HD179218 

2? 1.000 f^ 1 -^ ' ' ' ' 1 ' ' ' 1 ' ' ' ' ' ' ^ 

c - \ 

c ' \ 

e 0.100 r \ 

° ^^^"\ 

o V X 

D - S Ox 

c 

o 0.010 r \ - 

<D % 
N \ 

D \- - 

E - \| 
z 0.001 

0.6 1.2 1.8 2.4 3.0 

arcsec 

Fig. 2. Same as Fig. ^ to make a comparison with 
HD179218 for which there is no extended emission detected. 



extensions seen in the mid-IR range assumes a priori a 
disk-like geometry. 

Disks inclinations can be relatively easily determined if 
their emission are spatially resolved. However, HAe disks 
are generally dominated in the mid-IR by the innermost 
regions (1-30 AU). Our goal was to detect the emission 
from the intermediate regions of the disks (30-200 AU). 
The disk emission can be decomposed, in our image, in 
a central unresolved component plus an extended one 
(which geometry should reflect the true disk geometry 
at the distance scales achievable for our data). We first 
removed the "point-like" central emission component by 
subtracting a scaled PSF to the image of the object. 
The parameters of the point-like component (intensity, 
position) were computed automatically using a penalty 
functional (in order to avoid any visual bias) and then 
cross-checked visually using a dedicated graphical interface 
built in IDL. The resulting image, called residuals, is free 
from the central emission, so that the extended emission is 
enhanced and it is easier to determine the disk geometry. 
For each target, this processing was done using all available 
PSF measurements. The errors on the putative extensions 
were assessed when applying the same procedure to two 
PSF reference stars and, when possible, when comparing 
the extensions obtained from different nights. The results 
are shown in Fig. and |H1 

The first result is the elliptical shape of the extensions, 
which is characteristic of the emission of disks inclined with 
respect to the line of sight. An ellipse fit of the residuals 
gives an estimate of the disks inclinations and position an- 
gles. The results in terms of inclination and position angles 
are shown in Tab. [21 

The inclination of 33°± 5 found here for CQTau is in the 
middle ra nge of values found in the li terature, wh ich range 
from 63° ijTesti et al.ll2001|) to 14° (pent et al.ll2005[); in 
betw een we can find 48° ±£<, {PA = 105° ± 5) l|Eisner et alJ 
2004). 

For HD 163296, the disk inclination of 60°± 5° is in 
good agreement with that found by Mannings & Sargent 
(1997) of 58°. Conce r ning the position angle of the disk, 
iMannings fc Sartrentl (119971) fo und 126°± 3 with CO ob- 
servations, iGradv et al.l l)2000[) found 140°± 5 thanks to 
optical coronographic images. We found a value of 105°± 
10 when fitting our mid-IR data. This difference in the posi- 
tion angle could be related to the fact that our data are only 
sensitive to warm dust of which geometry could slightly dif- 
fer from that seen at shorter or longer wavelengths, if the 
disk contains for instance non axisymmetric structures. 



4. Inclination of the disk 

Extensions around two out of the three objects, namely 
CQ Tau and HD 163296. were already observed at other 
wavelengths. iTesti et "al"l l)2001l) have resolved the emission 
around CQ Tau at 7 mm and con cluded that i t was 
compatible with a disk-like geometry. IGradv et alJ (|2001f) 
have obtained coronographic images of HD 163296 with 
the Space Telescope Imaging Spectrograph on board the 
Hubble Space Telescope which revealed a circumstellar disk 
with a radius of 450 AU. Therefore, our modeling of the 



5. Modeling 

We have used a relatively simple parameterized model to 
investigate the dependence of the emission in the mid-IR 
on each parameter. We consider disks heated by irradia- 
tion from the central star. The density profile of the gas 
is parameterized as a function of r (radius) and z (vertical 
height above the disk mid-plane): 
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Table 2. Parameters of the elliptical fit of the image. In 
the first column, the inclination of the disk (0 degree means 
a faced-on disk). The position angle (PA) of the disk on 
the sky is measured counter-clockwise from the North. The 
distance corresponds to the maximum distance at which the 
disk is detected in the 20 /im observations. The brightness 
level in the last column corresponds to the brightness of the 
maximum distance at which the disk is detected here. 





i 


PA 


distance 


Brightness level 




(deg) 


(deg) 


(AU) 


(mJy/"2) 


HD 135344 


46 ± 5 


100 ± 10 


210 


170 


HD 163296 


60 ± 5 


105 ± 10 


215 


135 


CQ Tau 


33 ± 5 


120 ± 10 


290 


129 



and it is assumed that the dust is well mixed with the 
gas. The surface density is assumed to follow a power-law 
in radius: £(r) = £o(r/ro)~ p , with vq a fiducial radius. 
The scale height of the disk H p (r) is also assumed to be a 
power-law: H p (r) = Ho(r/ro) q . The inner radius of the disk 
(r=Ri n ) is located at the dust evaporation radius (1400 - 
1500K for silicate dust). The inner boundary (rim) is di- 
rectly exposed to the stellar flux and is puffed up since 
it is hotter than the rest of the disk. Here, we mimic the 
puffing-up of the rim predicted by DDN01 by a specified 
value of -Hp, i rl at r — Ri n , which is a parameter of the 
model. It should be noted that whether such an inner rim 
is indeed puffed-up i s still a matter of debate. Moreover, 
llsella fe Nattal i|2005l) have shown that the rim is proba- 
bly rounded-off due to the density-dependence of the dust 
sublimation temperature. This effect is not in cluded here. 

Fo r the dust opacitie s, we use those of iDraine ITU 
l|2001() . We use a MRN (jMathis et alJll977|) distribution 
of grains (n(a) cx a~ 3,5 ) with a size between 0.01 and 0.3 
/im. It is the disk surface layer which dominates the SED 
in the mid-IR range. The emission of the surface is made 
by small grains which trace the disk geometry. We will fo- 
cus in this paper on this component. Since the objects have 
no PAH (Polycyclic Aromatic Hydrocarbon) emission (or 
weak concerning HD 135344), we do not take into account 
this population of grains. 

Once the density profile is set, the dust opac- 
ities and stellar parameter s given, the code R ADMC 
iDullemond fe Dominild 120041 iPontoppidan fe Dullemondl 
2005) solves the temperature structure of the disk in 
a Monte-Carlo way using a variant of the algorithm of 
Bjorkman & Wood (1997). This Monte-Carlo code also 
produces the source terms for scattering, in the isotropic- 
scattering approximation. With a r ay-tracing tool (which 
is par t of the code RADICAL, see iDullemond fe Turollal 
(2000) for a detailed description) the SED and images 
can then be produced and compared to the observations. 
Comparative images are obtained by first resampling the 
maps to CAMIRAS sampling and then by convolving them 
with the PSF. 

6. Comparison model versus observation 

For each object, the best model shall fit simultaneously 
the SED (Fig. and the extension found at 20.5 (urn 
with our observations (Fig. EJ. Concerning HD 135344 and 
HD 163296, we used mainly the ISO spectrum to constrain 
the SED; as far as CQ Tau is concerned, IRAS photome- 



try and BASS points obtained bv iGradv et al J l)2005|) are 
used. Figure El shows the structure of the disk in terms of 
pressure and surface scale height. 

Multiple runs of the model are performed until a sat- 
isfactory fit to the observed spectrum and the extension 
found at 20.5 /jm is obtained. In the fitting procedure, 
the stellar parameters (T e ff, M*, R+, see Tab. |3J), the 
dust evaporation temperature (i.e. the position of the in- 
ner rim) at 1400 K, the outer radius and the inclination 
of the disk, the dust composition and size distribution are 
fixed. Other parameters, such as the pressure scale height 
for the inner rim (eventually puffed-up) and the outer pres- 
sure scale height (Ho), the power-law index of the pres- 
sure scale height (i.e. q, which has been fixed in the case 
of flared dis k to 9/7, a value determina ted by hydrostatic 
equilibrium IChiang fc Coldreichl the mass of the 

disk and the power-law index of the surface density (i.e. p) 
are estimated. Fitting the SED gives one solution among 
several degenerate combinations. A minimum value of the 
outer radius is derived from 20 /xm observations, which is a 
very strong constrain on the true disk size, thus removing 
largely the degeneracy on the set of model parameters. 

As first guesse s, we used those para meters found when 
fitting the SED bv lDominik et alJll2003ft for HD 135344 and 
HD 163296 and bv lChiang et al] l|2001j) for CQ Tau. Trying 
to fit the SED, we focused mainly on the Near-IR and mid- 
IR regions. Indeed, those regions are the regions where most 
of the reprocessed stellar energy (by the disk's surface layer) 
emerged, and therefore the most strongly affected by the 
model and the geometry of the disk. 

Not all the parameters are sensitive to the spatial 
distribution at 20.5 /<m (Tab.^J. The strong constrain that 
put the CAMIRAS images is the minimum outer radius, 
and we have to find a solution with all the free parameters 
to reproduce the shape of the SED, the extension and the 
total flux at 20.5 fim. 

For HD135344, the disk parameters deduced from fitting 
only the SED of HD135344 were close to those which allows 
to fit both the SED and our observations. 

For CQTAU, in addition to previous studies fitting only 
the SED, we had to take into account that the disk is quite 
extended at 20 fim (about 300 AU) and is observed with 
an inclination of 33° ± 5. The fast increase of the disk emis- 
sion in the [10-30] /zm range can be modeled only with a 
flared disk (with H° ut /R ou t > 0.1). In the framework of our 
modeling, we found that the only manner to obtain simulta- 
neously a quite low far-IR excess (as seen in the spectrum) 
with the observed extended 20 /xm emission is a low-mass 
disk of only 0.005 M Q with a pressure scale height of 58 AU 
at 450 AU. The total flux in the infrared excess compared 
to the stellar flux is determinated by the fraction of the 
central star energy intercepted by the disk. This covering 
fraction is linked to H° ut /R out , where the disk geometry 
thickness is maximum. Here, this parameter is the same as 
the previous study (energy conservation) but for a different 
outer radius. That means that the disk is less flared than 
deduced earlier. We obtain f or CQ Tau a disk mass 10 times 
small er than already found ijTesti et alJl2003t IChiang et alJ 
1200 1^) and this mass only traces the small grains. The disk 
emission at mm wavelength is determined by the mid-plane 
grain and disk properties, and it is not affected by the na- 
ture of the surface dust. Consequently, in this paper, we 
have not tried to select the best parameters for the mid- 
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plane dust to fit mm observations. Our underestimated flux 
at mm wavelength suggests that there must be big grains 
in the disk mid-plane of CQ Tau to recove r the measured 
flux at these wavelengths l|Testi et al.ll2003fl . 

HD i63296 is classified as a group II object (interpreted 
as being surrounded by flat or a self-shadowed disk); it 
should be in principle much more difficult to resolve the 
disk in the mid-IR range. Surprisingly, some extended mid- 
IR emission at 20 /im, although less prominent than for 
group I objects, is however observed. In th e modeling, we 
used a s first guess the parameter found in iDominik et al.l 
( 2003|, who modelized the SED with a flared disk cut at 50 
AU. Here, we modified the large-scale pressure scale-height 
in order to mimic a weakly flared disk. The SED is well re- 
produced with a disk having little flaring (Fig0]and Tab.[3J) 
while the extension at 20 fim (Fig. |3J constrains the disk to 
have a minimum outer radius of 200 AU with a disk mass 
of 0.01 M . 

7. Conclusions and future work 

We have shown the strength of mid-IR imaging to constrain 
the disk properties. There are still a very limited number 
of objects with extended emission spatially resolved and 
it is not yet possible to draw statistical conclusions about 
the spatial structure of disk around HAe stars. With the 
advance of mid-IR i nstruments on 8 me ter class telescope, 
such as the VISIR (jLaeaee et al.ll20?)^l instrument avail- 
able on the MELIPAL Very Large Telescope (VLT) at the 
European Southern Observatory (ESO), higher angular 
resolution will be available and the field will further 
develop; more quantitative studies will be possible and 
the goal of retrieving detailed disk surface density profiles 
from the observations should be achievable. 
Observations should not be limited to the 20 /im atmo- 
spheric window; observations in the 10 /zm atmospheric 
window are also very promising, especially for those HAe 
stars whose spectrum features the so-called PAH bands 
features at 3.3, 6.2, 7.7, 8.6, 11.2, 12.7 ^m. PAH bands 
are attributed to vibrational relaxation of UV-pumped 
Polycyclic Aromatic Hydrocar bon molecules containing 
about 50-100 carbo n atoms ijAllamandola et ahl [l989; 
iPueet fc Legerl Il989|) . Their emission, as a function of 



HD1 35344 mode 



the distance r to the star, drops with a 



power-law, 



much slower than thermal emission from large grains in 
thermal equilibrium. PAHs emission is thus a promising 
probe to study flaring disks at large distances from the 
star, with the good angular resolution achieved now on 
large ground-based telescopes at 10 microns. Fitting by 
sophisticated models a combination of interferometric 
observations in the near- and mid-IR, which probe the 
inner-most disk regions, with single dish observations in the 
near and mid-IR which probe intermediate disk regions, 
is clearly the way to clear up the field in the next few years. 



JJ 0.010 




1 . 8 

arcsec 



Fig. 3. Annular average intensity aperture normalized to 
the peak for the observations (full line), model (dashed 
line) and the PSF associated (dot-dashed line). We have 
constructed an image with the DD04 model that we have 
convolved with the PSF of the night associated. The error 
in the profile comes from the noise RMS in the image. 
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Table 3. Stellar properties and fit parameters. The stellar p arameters of CQ Ta u ar e taken in | Chia ng et al.] 1I2OOI). For 
HD 135344 and HD 163296, the parameters are taken from lMeeus et all (|2001fl and Ivan Boekel et al l l)2005|) . H™/R in 
characterizes the puffed inner rim. q is the power-law index of the scale height and p that of the surface density. H°"*/i? out 
characterizes the flaring angle of the disk. The last column is the flux at 20.5 ^tm calculated from the modeled image. 



Object 


distance 


M* 


Teff 


R* 




Hp™/-Rin 


Rout 


H°"*/ Rout 


P 


q 


i 


F(20.5 ^m) 




(PC) 


(M ) 


(K) 


(Rq) 


(AU) 




(AU) 








(degrees) 


(Jy) 


HD 135344 


84 


1.3 


6750 


2.1 


0.24 


0.065 


800 a 


0.21 


0.8 


9/7 


60 




(DDN03) 
























HD 135344 


140 


1.3 


6750 


2.1 


0.24 


0.065 


200 


0.12 


0.8 


9/7 


45 


5.7 


(in this study) 
























HD 163296 


122 


2.5 


10500 


1.7 


0.45 


0.033 


50 


0.07 


0.2 


9/7 


65 




(DDN03) 
























HD 163296 


100 


2.5 


10500 


1.7 


0.45 


0.033 


200 


0.05 


1. 


1.07 


60 


16.6 


(in this study) 


























CQ Tau 


100 


1.7 


7130 


1.3 


0.23 


0.0002 


180 


0.13 


1.5 


9/7 






(Chiang2001) 
























CQ Tau 


100 


1.7 


7130 


1.9 


0.23 


0.018 


450 


0.13 


0.3 


9/7 


33 


20 


(in this study) 

























a Disk size compatible with new measurement due to the re-evaluation of the distance. 



Table 4. Effects of the free parameters of the model on the SED and spatial distribution of the 20.5 /j,m emission in the 
CAMIRAS image. 



Parameter 


Influence on SED 


Spatial distribution at 20.5 /im 


So 


Influence on the total flux in the mid-IR 
and more strongly for A > 60 /im 


No influence 


P 


No influence in the mid-IR 


little influence 


Hp n /-Rin 


modify contrast mid/near-IR 
(structure inner rim/shadow) 


No influence 

(inner rim + shadow in the first pixel) 


1 


modify the whole shape of the SED 


large influence 


Hp"'/ Rout 


modifiy flux in the mid and far-IR 


little influence 


Rout 


modify far-IR emission 


large influence 
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Fig. 4. Modeled spectrum (dashed line) between 1 and 100 
Urn. For HD 135344 and HD 163296, the ISO spectrum is 
overploted in bold line. For CQ Tau, there are some photo 
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Fig. 5. Modelized vertical structure versus the radius. The 
pressure height scale (H p /R) is in full line and the surface 
metric" points "fror dGradv" et"aTI ijfMT we have" added" "the nei S ht scale ( Le - the hei S ht of the disk photosphere above 
one from this study at 20.5 /an and IRAS points at 60 and the mid-plane H s /R) is the dashed line. 
100 ^m (full circle). 



10 



C. Doucet et al.: Mid-InfraRed imaging of Herbig Ae stars 



HD1 63296 

0.193 Jy 




-0,045 



Fig. 6. HD 163296 extended emission (central point source 
removed). West is up, North on the left. The pixel size 
is 0.29 arcsec. The disk has a surface brightness of 0.59 
Jy/" 2 (S/N=27 for a ~ 0.002 Jy) for the brightness part 
in the direction east/west and 0.23 Jy/" 2 (S/N=ll) for the 
less bright part in the direction north/south. The noise is 
calculated using a sigma-clipping technique. 



CQ Tau 

0.082 Jy 




-0.019 



Fig. 7. CQ Tau extended emission (central point source 
removed). The same orientation as FigH3 The pixel size is 
0.29 arcsec. At 0.9" from the center, the surface brightness 
is 0.94 Jy/" 2 (S/N=15 for a ~ 0.004 Jy). At 2" from the 
center, the surface brightness is 0.23 Jy/" 2 (S/N=5). 



HD1 35344 



0.021 Jy 




-0.017 



Fig. 8. HD 135344 extended emission (central point source 
removed). The same orientation as Fig EI The pixel size is 
0.29 arcsec. At 0.9" from the center, the surface brightness 
is 0.25 Jy/" 2 (S/N=15 for a ~ 0.0015 Jy). 



